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THE DEFICIENT QUARTIC SPLINE INTERPOLATION 
Y.P. DUBEY, R.K. DUBEY AND S.S. RANA 


Abstract: In the present paper, we have studied the existence, uniqueness and convergence properties of 
deficient quartic spline interpolation. 


1. INTRODUCTION 


Quartic and higher degree splines are still popular for best approximation (see 
de'Boor [1]). In the study of lower degree piecewise spline functions there is a 
disadvantage, that we have corners at joint of two lower degree pieces and therefore to 
achieve a prescribed accuracy more data than higher methods are needed. Rana and 
Dubey [9] generalized the result of Garry and Howell [7] and obtained best error 
bounds for quartic spline interpolation. Dubey and Shukla [4] have obtained error 
bound for quartic spline interpolation. For various aspects of cubic, quartic and spline 
of degree six reference may be given to Meir and Sharma [8], Hall and Meyer [6], 
Gemling-Meyling [5], Dubey [3]. In this paper, we have investigated the existence 
uniqeueness and convergence properties of deficient quartic spline interpolation which 
match the given function at values at mesh points, interior points with appropriate 
boundary conditions. 


2. EXISTENCE AND UNIQUENESS 


Consider a mesh P of [0,1] given by 0=x< x, < sm <x, =1 such that 


x- x, =h,i=0,1,...n-1 and h=max h,. Let s(4,P) denote the set of real algebraic 


i+l 


polynomials of degree not greater than 4 i.e. 
S(4,P)={s,,s,eC [0,1],i=1,2...n-1}. 


where s,is the restriction of s over [x,,x,,,]. Let S* (4,P) denotes the class of all 


i+] 


deficient quartic splines S(4,P) which satisfy the boundary conditions, 


s" x)= f"); 
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s"(x )=f"(x,). (2.1) 
We shall prove the following. 


THEOREM 2.1 : Let f', f” exist. Then there exist a unique deficient quartic spline in 
S* (4, P) which satisfies the interpolatory conditions, 


s(x,)= f (x,) i=0, 1, 2......n (2.2) 
s(a@,)= f(a@,) i=0, 1...., n-1 (2.3) 
where a,=x,+Oh, if o< <3-Ẹv6 and < nh, > 


is non increasing. 

PROOF : Consider a quartic polynomial P(z) on [0,1]. We can easily verify that 
P(z)=P(0)q,(z)+ P(1)q,(z)+ P(0)q, (2z) 

+ P"(O)q,(z)+ P"A)q,; (z), (2.4) 
where 

q,(z)=14+ KC -60 -6+120°)z-12 (0 -1)z°+6(0 -1)z*]/64A, 

q,(z)=[(60* -12 0°)z+12 62° - 60z*]/6A, 

q,(z)=[6z -12 z°+6z*]/64A, 

q, (z)=[(20f+50°-30°)z+3(0* +0 -20°)z’” 


+(-0-50+60°) z+(0° -30° +20)z*]/6A, 





q,(z)=K-0*+0°)z+(0*-0)z°`-(0°-0)z*]/64A, 
with A =(0°+60-20°) 


- x,) 





Now writing += 2 , (2.4) may be expressed in term of restriction s, of s as 


h, 


follows :- 


s (x)= f(x), + fx) A+ f(@)a,0+¢,05,,4,)+4,@)5,(x) (2.5) 
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which clearly satisfies the condition (2.1) - (2.3) and s(x) is a quartic in [x,,x,,,] for 
i=0,....n. Since sec’ [0,1], so applying continuity condition of first derivative of s in 
(2.5) we have 


3 


30° -@)s,_,+{h,_,(20° -30° +80) 





+ h, (30° -20°-50°)} 5, +h, (0°-0°)s,,,=F, (Say) 


1 





T = ae 3 a 2 
where Pic 60 -6+120) 0(0 -20 P 


ði 








LL h, hy, J 
(-60°+6+120°-120)f,, 6 
f&a) 
hi h 
6 60° (0-2) ] 
ee f Gag) (2.6) 
i | 


Clearly the coefficient matrix of the system of equation (2.6) diagonally 
dominant under the conditions of Theorem 2.1 and hence invertible. This complete 
proof of theorem 2.1. 


3. ERROR BOUNDS 


Following the method of Hall and Meyer [6] in this solution, we shall obtain 
bounds of error function e(x)= f(x)- s(x) for the spline interpolant of Theorem 2.1 


which are best possible. Let s(x) be the second time continuously differentiable quartic 
spline function satisfying the conditions of theorem 2.1. Now considering 
fec{o,1] and writing L,[f,x] for the unique quartic which agree with 


f(x,) f(x,,,) f(@,) f"(x,)and f"(x,,,), we see that for xe[x,.x,,,],we have 





| FO) = s(x) [=] f= 5,0) | 
<| FOA Ef xt | LIF, x1- 5,00] (3.1) 


First we have to obtain the bounds of right hand side of (3.1) by Cauchy Theorem 
(See Davis [2]). Thus, we have 


5 


h; 
| f()LIf,xI|< a t) (t-9)| (3.2) 
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jiu 
where r= 





Xx. 
~ and F =max | f°’(x)|- 
O<x<l 


To get the bounds of | L,[ f, x]- s,(x) |,we have from (2.4) 


Alf. xl- s, (x) =h ex), (+h e" (xq, AN (3.3) 
Thus, 


[|LUf.x]-s i{(x)| 


<| hy e"(x,,) 4, +h, e"(x,) 4.1) | (3.4) 


Let max | e(x,)| exists for i= j then the inequality (3.4). 


[LL f.x1-s,(x)| 
<h le" lla, O+ ON 


Let k(t)=max [0|0°+60-3|+3|0°-26°4+1| t+6|8 -1| +30-0)t°]. 


0<0<1 


Thus | L [f,x]-s;(x)|< h‘ |e" (x,) |k (t) (3.6) 


Now, we proceed to obtain |e"(x,)| replacing s"(x,) by e"(x,) in equation (2.6), 


we have 


h, (30° -30° +30 -@)e, 





{h (20° -30° +0) 

+h,(30° -20° -50°)} e; +h (80 -1)0° e, 
4 3 2: S 

=F, -h,, (50 -30°+30°-6)f, , 


-{h,,(20° -30° +0)+h, (30° - 20° -50°)} 


f, —h, (0 -1)0° f,,=E(f) (say). (3.7) 


In view of that E(J is linear functional which is zero for polynomials of degree 4 
or less we can apply Peano Theorem (see Davis [2]) to obtain 


5) 


Dales o (3.8) 





( 
E= $ 


¥j-1 4! 
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Thus, from (3.8) we have 


1 Xj 
CES T Eaa (3.9) 


Where F =mx | f *(y)]. 


Further, it can observed that from (3.9) in x,,<x<x,,,we get 








F -60 -6+120° 60(0° -280° +12) i 
E{(x,- y), J=l{ }(x,-y), 
h, h, 
6 6 60° (0-2) 
m a Ree sae aa P 
h, h, h, 
-12 {h,_,}20* -30° +0) 
+h, (30° -20 -50°)(x,- y)? -12 h, (0f -@°)(x,,,-y).+ (3.10) 


In order to estimate the internal of r.h.s. of (3.10) we rewrite the above 


expression in the following symmetric form about x, to see 





60° 0 -2Xx,-y+h,) =2h (0-1) [x,-y+h,]? AS Hea, 


6 
=—[(9* -20° +1)(x,- y) +(40*-80°+0)h,(x,- y)+20° 
h, 


i 


IA 
IA 
Q 


(20° -50 +3)h?] Kx,- y)’] xS ysa, 


6 4 3 2 3 2 
= —[(@ -20 +1)(x,- y) +0(40 -80 +1) 
h 


i-l 


h, (x, — y)+20° (20° -5043)h,,(x,-y)’, dS pee 


i-l i 





60 2 2 2 
m (C8 —2)(x,-y-h,_,) +2h (@-D}lx,-y-A_ J, 


Thus, it is clear from above expression that £[{(x- y)] is non negative for 


x,,<x<x,,, given by 


i~} 
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fl E(x- y)* |dy = (nf +h?) (3.11) 


3 


[o | 
where EJE A -124 6° + 400 0° — 543 0° ab 
2 


Combining (3.9) and (3.11), we have 


n°). (3.12) 


i-l 


K (0); 4 
|E, Is F lh, 


Thus, from (3.8) and (3.11) 











, , FK (0) (hf, +h-) 
le (x,)|Sle, |< ; 3 : : 7 : ; (3.13) 
51[(-30° -30° + 20)h_, +(60°+0* -30° )h,] 
Now using (3.2), (3.7) along with (3.13) in (3.1) we have 
he 2 2 2 Pt 
Je(x) |< — |t’ A-t)? (t-0)|F +h le (x) |k) F 
5! 
h° FK (0)h°k 
=—|1°(-1)*(t- 0) |F + — me = (3.14) 
5 (60° — 20° -60° +20) 
h° 
<—F|C(t)| (3.15) 
5! 


where c (1)=|1?—1)°(@-0)+K. (@)K(1)| 


K (8) 
(60° -20 -60° +20) 





and K, (0)= 


Thus we prove the following. 


Theorem 3.1. Suppose s(x) is the quartic spline of Theorem 2.1 interpolating a 
function f(x) and fec’[0,1]. Then 


5 


h 
Lore IW], (3.17) 


where K = max |C (t) |=max |t (1—-1°)(t-0)+K, (@)k(t) |, 


O<r<l O<t 


and 


7/10 
THE DEFICIENT QUARTIC SPLINE INTERPOLATION 


5: 
|e"(x,)|s 





K | (0) max IF]. (3.18) 


5! 


Further more, it can be seen easily that K in (3.17) can not be improved for an 
equally spaced mesh. Inequality (3.18) is also best possible. Also, we have 


5 


h 
Je" (x, ) |SK, — mx | fW] , (3.19) 
5! O<x<l 


where «xK is positive constant. 


Equation (3.15) proves inequality (3.17) where as inequality (3.18) is direct 
consequence of (3.13). 


Now we turn to prove that the inequality (3.17) is best possible in the limit. 


5 


Considering f (x)= * and using Cauchy formula Davis [2] we have 
5! 


T 5 il >: h? 2 1- 2 -9 
A E E eee (3.20) 
‘| 5 | 5! 5! 





For the function with equally spaced knots from equation (3.7), the following 
relation holds 





Px? | 3 

se es -(50°-30°+30°-@)e , 

| st | 5! l 

+(0 -80° +367 Je +(@-1)0°e",, (3.21) 


Consider for a moment 


K (O)h° 





Cu 
J 


mes ane (3.22) 
60° 0-0) ` l 


We have from (3.4) 


K (0)h° 





s(x) -L,[f, x] (¢,(t) +q; (t)) 


610? aoe 


K, (@)h° 
ME gee? +60 —3)t+20(20°+1-36")t° 
610° (1-6)° 


+60(0 -1)t +30(1-6)t'}. (3.23) 
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Now combining (3.20) and (3.23) we get for 


K (0 
ues {0°(0 +60 -3)t 





ie 
s(x) = f (x) =—I[— - 
5! 60 d-0) 


+20 (20° +1-30°)t? +60(6-1)t° 
+301 -0)t'}+t A-4 0-0). (3.24) 


From (3.24) itis clearly observed that (3.17) is best possible, provided we would 
prove that 


k (O)h° 
6!10 0-0)? 





(3.25) 


i-l 1 Cia 


In fact (3.25) is attained only in the limit, the difficulty will appear in the use of 
boundary condition ie. e"(x,)=e"(x,)=0. 


However, it can be shown that as we move many sub-intervals away from the 
boundaries 


K (@)h° 


e" (x,) > ——— 
6!10 (1-0) 


. For that we shall apply (3.20) indutively to move away 


from the end conditions e"(x,)=e"(x,)=0. 


The first step in this direction is to establish that e"(x,)=0 for some i=1,2,..n 


which can be shown by contradictory result. Let e"(x,)<0 fori=1,2,...n-1. 


Now making use of (3.18) we have 


K (0)h° 


2 





1 
> max |e"(x,)|>—[(50' -30° +30 
2 


2 2 O)e,, + (0 - 1)0°e,, 
6'0 (1-6) OSes 


1 


1 4 3 2 n 3 2 p Bro 
>—[(50 -30 +30 —@)e,,+(0—-80 +30 Je, +(0-1)0 e, 
2 


1 


11K, (ayn? | 
|e | 


2j 5! | 


This is contradiction, hence e"(x,)>0 


limiting case e"(x,)> 


K (0)h° 
60(1-0) 
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Now from (3.21) 


K (O@)h- 





(0 -80°+30)e = a (50° -30°+30° -@) 
e,,-(@-1)0°e,,, 
E Oe 
5! 
1 K (8)h° 


e, < 
(0-80°+30) 5! 





Now again using (3.26) in (3.21) we have 


3 
K (0)h 
(0 -80° +30°)5 





e"(x,)< 





2 
[  60t-40°-0 (66'-40°436°-6 | 
<l1 ; >t ; ; + 
| (@-80° +30 0 -80° +30 | 


Now it can be seen easily that r.h.s. of (3.28), > 


3 


K (O)h° 
60(1-0) 


K (O)h° 
601-0)? 


which verify proof of inequality (3.19). 
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(3.26) 


(3.28) 


and hence in the 


Thus, corresponding to the function fœ)=> (3.27) and (3.28) tends to 
5! 


~in the limit for equally spaced knots. This complete proof of Theorem 3.1. 
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